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\£) , Abstract. The supernova-driven interstellar medium in star-forming galaxies has Reynolds 
numbers of the order of 10 6 or even larger. We study, by means of adaptive mesh refinement 
hydro- and magnetohydrodynamical simulations that cover the full available range (from 10 kpc 
to sub-parsec) scales, the statistical properties of the turbulent interstellar gas and the dimension 
of the most dissipative structures. The scalings of the structure functions are consistent with a 
log-Poisson statistics of supersonic turbulence where energy is dissipated mainly through shocks. 
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. 1. Introduction 

The interstellar medium is a highly turbulent compressible flow with Reynolds num- 
bers of the order of 10 6 or even larger (e.g., Elmegreen & Scalo 2005). Supernovae are the 
main source of energy of the interstellar gas, which is dominated by nonlinear processes, 
including heating and cooling, that transfer energy among a wide range of scales. It is un- 
clear what are the driving scales in the interstellar medium as observations show different 
values ranging from a few parsec to tens or even hundreds of parsec (see, e.g., Haverkorn 
et al. 2005). Furthermore, it is still under debate which dissipation process dominates 
the interstellar turbulence. High-resolution simulations of the interstellar medium that 
cover a wide range of scales, including the injection and viscous scales, are valuable tools 
to adress these issues. We, therefore, carried out kpc-scale high resolution (up to 0.625 
pc) hydro- and magnetohydrodynamical simulations of the supernova-driven interstellar 
medium (Avillez & Breitschwerdt 2004, 2005, 2006) allowing us to tackle simultaneously 
the disk and halo evolution and to understand the issues referred above. In this paper 
we briefly discuss these results. Section 2 deals with the injection scales, the scalings of 
the velocity structure functions and dissipation processes are referred to in section 3, 
followed by a discussion of the results and their implications (section 4). 



2. The Injection Scale of Interstellar Turbulence. 

The outer scale of the turbulent flow in the ISM is related to the scale at which the 
energy in blast waves is transferred to the interstellar gas. Such a scale can be determined 
by using the so-called two-point correlation function R^j [l, t) — (iii(x+l, t)uj(x, t)), where 
Ui are the components of the fluctuating velocity field u. The diagonal components of Rij 
are even functions of / and can be written in terms of dimensionless scalar functions f{l, t) 
and g(l,t), which satisfy /(0) = g(0) = 1 and /, g ^ 1, that is, Rn/u 2 = f(l,t) and 
Rij/u 2 = g(l, t) if i = j 7^ 1 and zero if i ^ j with u = (i(w)) 1 / 2 . In these equations R\\ 
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and i?22 = R33 are known as the longitudinal and transverse autocorrelation functions, 
respectively. 
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Figure 1. History of the characteristic size (given by Lii) of the larger eddies (left panel) and 
of the ratio L22/L11 (right panel) for the HD (dashed lines) and MHD (solid line) runs for 1.25 
pc (black) and 0.625 pc (red) resolutions. 

The outer scale of the flow is given by Lu — J + °° f(l,t)dl, which is calculated in a 
region with a linear size of 500 pc at a distance of 250 pc from the edges of the compu- 
tational domain such that we avoid the periodicity effects of the boundary conditions. 
The top panel of Figure [1] shows the history of Lu in the last 50 Myr of evolution of the 
simulated interstellar medium for 1.25 (in black) and 0.625 pc (in red) resolutions. The 
average integral length scale in the three simulations is 73-75 pc and seems independent 
of the resolution. In all the cases there is a large scatter of Lu around its mean as a result 
of oscillations in the local Galactic star formation rate, where the formation and merging 
of superbubbles, is responsible for the peaks observed in the two plots. The correlation 
length scale measured in the simulations is remarkably similar to that determined by 
Kaplan (1958) who found that interstellar velocities start to decorrelate at a scale of 80 

PC 

The transverse integral length scale given by L22 — J Q g(l,t)dr equals 0.5Ln in the 
case of isotropic turbulence. In the present simulations we have 0.2 < L22/L11 < 1-3 
(bottom panel of Figure [T]). In spite of the large scatter the time average of the L22/L11 
over a 50 Myr period is 0.51 and 0.6 for the HD and MHD runs, respectively. The 
discrepancy of < L22/L11 >~ 0.5 by about 20% in the MHD run is a consequence of 
the anisotropy in the flow introduced by the magnetic field. Overall the field strength is 
still too low to produce a larger deviation. In case of the HD run < L22/L11 >~ 0.5, 
indicates that in a statistical sense the interstellar unmagnctiscd turbulence is roughly 
isotropic. 




3. Velocity Structure Functions and the Most Dissipative Structures 

Turbulent flows are usually described in a statistical fashion by the velocity structure 
functions of order p defined as S p (l) — {Svf) oc l^ p >, where ((p) is a scalar, 5vi — 
\v(x + I) — v(x)\ with v(x + I) and v(x) being the velocities along the a;— axis at two 
points separated by a distance / such 77 <C I <C Lu, with Lu being the energy integral 
scale and rj = v 3 / 4 e~ 1 / 4 the Kolmogorov microscale, respectively. Here, v is the kinematic 
viscosity and ( ) stands for the ensemble average over the probability density function of 
Svi. From the Kolmogorov (1941; hereafter denoted by K41) theory one derives 

(Svf) oc {8v?}^/« 3 \ (3.1) 
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with I 3> rj. This property, referred to as extended self-similarity (ESS), indicates that 
S p is dictated solely by the third-order structure function S3 via the set of scalings 
C(p)/C(3) and is valid in a wide range of length scales for large as well as small Reynolds 
numbers even if no inertial range is established (Benzi et al. 1993). Using the ESS concept 




Figure 2. Left: Exponent C(p)/C(3) for the structure function versus order p. The full line 
corresponds to £(p) = p/3 (K41 theory), bullets correspond to data of Benzi et al. (1993), red 
and blue lines refer to the She-Leveque (1994) and Burgers-Kolmogorov (Boldyrev 2002) models, 
respectively. Triangles represent data from HD simulations with 1.25 pc (black) and 0.625 pc 
(red) resolutions and the green squares refer to the MHD run with 1.25 pc resolution. Right: 
Fractal dimension D of the most dissipative structures (derived from eq. I3.2|l as function of the 
exponent C(p)/C(3), with p = 2,4, 10, shown in the left panel for the three simulations. 

we determined the velocity structure functions (Svf), with order p = 2, 4, • • • , 10 as a 
function of (Svf) and their best fits slopes C(p)/C(«0 ( see Avillez & Breitschwerdt 2006) 
for the simulated interstellar gas. The slopes are displayed as function of the order p in 
the left panel of Figure [21 where triangles (black corresponds to Ax = 1.25 pc and red 
to Ax — 0.625 pc) and squares refer to the HD and MHD runs, respectively. 

From these slopes one can determine the Hausdorff dimension of the most dissipative 
structures through the solution of (Dubrulle 1994) 

H = e(i-A)p + 3 A_(i-/? ep )» (3-2) 

where A = 1 — O is a parameter depending on the Hausdorff dimension D of the most 
dissipative structures, /3 = 1 — A/(3 — D) is a measure of the intermittency and is 
a parameter that depends on cascade of the energy transfer in the inertial range. If 
this cascade is Kolmogorov, then 9 = 1/3. The above mentioned ratios correspond to 
D = 1.9 — 2.02 (right panel of Figure implying that the energy injected by supernovae 
into interstellar turbulence is dissipated preferentially through shocks (i.e., 2D surfaces). 

Further analysis on the Hausdorff dimension of the most dissipative structures can 
be drawn from the comparison of the simulated and theoretical predictions (eq. 13. 2ft of 
the C(p)/C(3) variation with p shown in Figure [2] (left panel). For comparison we also 
show the variation of the ratio C(p)/C(3) observed experimentally by Benzi et al. (1993) 
and predicted by the K41 and She-Leveque (1994; hereafter denoted by SL94) models in 
incompressible turbulence and the Burgers-Kolmogorov model (Boldyrev 2002; hereafter 
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denoted by BK02) for supersonic turbulence. The K41 model considers a log-normal 
statistics for the transfer of energy from large to small scales and has no corrections to 
intcrmittency. The SL94 and BK02 models use a log-Poisson statistics (see discussion in 
Dubrulle 1994) to describe intermittency. 

The figure shows that there exists a small discrepancy between the MHD data (green 
squares) and the BK02 curve and a decrease of D with increase of p. This variation, 
although within the errors of numerical noise, indicates a tendency toward filamentary 
dissipative structures due to the anisotropy induced by the magnetic field. 

4. Discussion and Final Remarks 

The dissipation of energy in interstellar turbulence proceeds through shocks and the 
variation of C(?>)/C(3) with p is most consistent with a log-Poisson model where at the 
inertial range the cascading of energy behaves more like Kolmogorov. Similar scalings 
have been found by Boldyrev et al. (2002) and Padoan et al. (2004) for the dissipation of 
energy in molecular clouds. Note, however, that our simulations include molecular clouds 
and its fragmentation and cover a larger range of scales (from kpc to 0.6 pc) than those 
used by these authors. The simulations also show that only a small amount of energy 
is dissipated through shocks in the inertial range and the energy decay follows a quasi- 
Kolmogorov spectrum. However, the shock structures start to play an important role in 
energy transfer and dissipation near the small scales. 
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